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ABSTRACT 

We study the sources of biases and systematics in the derivation of galaxy properties 
of observational studies, focusing on stellar masses, star formation rates, gas/stellar 
metallicities, stellar ages and magnitudes/colors. We use hydrodynamical cosmolog¬ 
ical simulations of galaxy formation, for which the real quantities are known, and 
apply observational techniques to derive the observables. We also make an analysis 
of biases that are relevant for a proper comparison between simulations and observa¬ 
tions. For our study, we post-process the simulation outputs to calculate the galaxies’ 
spectral energy distributions (SEDs) using Stellar Population Synthesis models and 
also generating the fully-consistent far UV-submillimeter wavelength SEDs with the 
radiative transfer code sunrise. We compared the direct results of simulations with 
the observationally-derived quantities obtained in various ways, and found that sys¬ 
tematic differences in all studied galaxy properties appear, which are caused by: (1) 
purely observational biases, (2) the use of mass-weighted/luminosity-weighted quan¬ 
tities, with preferential sampling of more massive/luminous regions, (3) the different 
ways to construct the template of models when a fit to the spectra is performed, and 
(4) variations due to the use of different calibrations, most notably in the cases of 
the gas metallicities and star formation rates. Our results show that large differences 
can appear depending on the technique used to derive galaxy properties. Understand¬ 
ing these differences is of primary importance both for simulators, to allow a better 
judgement on similarities/differences with observations, and for observers, to allow a 
proper interpretation of the data. 

Key words: galaxies: formation - evolution - cosmology: theory - methods: SPH 
simulations - SPS models - radiative transfer 


1 INTRODUCTION 


In recent years, large galaxy surveys such a s the 2dFGR S 
(Two-degree-field Galaxy Redshift Survey, Colles^ 

SDSS (Sloan Digital Sky Survey, lAbaza-iian et al.l 1200^1 
and 2MASS (Two Micron All-Sky Survey, Skrutskie et al.l 
I2OO6I) . have opened up the possibility to statistically study 
the properties of galaxies in the Local Universe, reveal¬ 
ing their great diversity: even for a narrow range in stel¬ 
lar mass, galaxies appear in a large variety of morpholo¬ 
gies, gas fractions, star formation rates (SFRs) and chem¬ 
ical abundances. These observations have also allowed to 
identify important relations such as the ma ss-metallicity 
(iGarnbt fc ShieidJi 19871 : iTremonti et aLll2004l) . and to mea¬ 
sure the corresponding scatter which encodes relevant in¬ 
formation on the galaxies’ evolution. This wealth of data 
gives important insight on the process of galaxy forma¬ 
tion and evolution, revealing the action of physical mech¬ 


anism s occurring i n galaxies, both internal - e.g. fe ed- 
back (iFabiarJ I2OI2I ) , cooling llThoul fc WeinberJ Il 995 l) - 
and in relation to larger-scale mechanisms - m ergers 
llNaab et aJj 200jjNajAGc4ianssqri_j£U3strjjmi^^ inter 


actio ns (lScuddeji.et_alJ 2012: Stierwalt et al. 20151 ). accre¬ 
tion (IPutman, Peek fc Jounall2012l) . All these leave imprints 


on the shape of the spectral energy distributions (SEDs) 
which constitute the primary source of information from 
large galaxy surveys. In fact, in recent years it became possi¬ 
ble to obtain the full SEDs of galaxies at wavelengths from 
the X-Ray to the radio. In particular for galaxy studies, 
wavelengths from the ultraviolet to the far infrared are the 
most relevant as the light coming from the stars and inter¬ 
stellar gas/dust dominates the spectra in this wavelength 
range. 


In addition to observations, numerical simulations 
are also a useful tool to study the formation and evo- 
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lution of galaxies and to investigate the links between 
a galaxy’s formation, merger and accretion history and 
its final properties. Understanding these relations is rele¬ 
vant for a correct interpretation of observations of galax¬ 
ies at different cosmic epochs, and for a reconstruc¬ 
tion of the galaxy’s histories. Progress on the simulation 
of realistic galaxies has been presented in many recent 
works (e.g. Gnvernato_et_alJ 2007; ScannM)ieco_et_al, 200^ 


Aurner et alT 2013ir Voee lsberger et al.l l2014l : I Nelson et all 

2 OI 5 I : ^have et al.l l2015h that, together with advances in 
computational resources, are starting to allow the simula¬ 
tion of relatively large cosmological volumes, recreating a 
virtual universe where galaxies are naturally diverse as the 
result of their individual evolution. 


Despite the individual progress in observations and sim¬ 
ulations, much uncertainty remains in the comparison be¬ 
tween them even when these are used to decide on the 
successes and failures of the models and are recognized 
as an important aspect of the interpretation of observa¬ 
tional results. In fact, previous works using simulations 
have shown that, as the methods usually applied to de¬ 
rive the properties of simulated galaxies are very different 
from observational techniques, several biases mi ght be intro¬ 
duced making the comp a risons u nreliable (e.g. Abadi et al.l 
200^^ Gnvernntoetal _ 2009|;_ Scannapieco et al.l 201^ 


Snyder et al. 201ll : Munshi et ahl'20131: Chri stensen et al.l 

20141 . For example. IScannapieco et al.l (120101 1 showed that 
disc-to-total ratios can vary up to a factor of 250% when 
analysed using kinematic or photometric disc-bulge(-bar) 
decompositions. In addition, some authors pointed out the 
effects of the stochastic sampling of the Initial Mass Func¬ 
tion (IMF) in observations, in particular when properties are 
deri ved fitting Stellar Pop ulations Synthesis (SPS) models 
(e.g. iKrumholz et al.|l2015h . In order to properly judge the 
agreement between simulations and observations and make 
it possible to better understand observational results, it is of 
primary importance that these biases are well understood. 


In this work, we do a systematic analysis of the sources 
of biases in the comparison between simulated and observed 
galaxies, creating synthetic SEDs for our simulated galax¬ 
ies and then deriving the galaxy properties using differ¬ 
ent techniques that mimic those used in galaxy surveys. 
As the true values for properties such as SFRs, metallici- 
ties, ages and stellar masses are known in the simulations, 
we can quantify how close the values obtained observation- 
ally are to the real ones. To create the synthetic SEDs, we 
use three methods with increasing complexity: first, we use 
dust-free SPS models; second, we include dust using a sim¬ 
ple parametrization; and third, we calculate the transfer of 
the galaxy light through the Interstellar Medium (ISM). In 
this last case, we post-process the simulation outputs with 
the 3D p olychromatic Monte Carlo ra diative transfer code 
SUNRISE (Ijonsson. Groves fc Coxl|2009l l. which simulates the 
propagation of photons in a dusty ISM. Using the synthetic 
SEDs, we derive stellar masses, stellar ages, stellar and gas 
metallicities and SFRs. By comparing the direct results of 
the simulations with the observationally-obtained quanti¬ 
ties, we study (i) how reliable and meaningful simple com¬ 
parisons between simulated and observed galaxies are, {ii) 
how the different assumptions involved in the creation of the 
SED affect the derived galaxy properties, and (Hi) how we 


can reliably test the agreement between observations and 
models. 

In this paper, we describe the techniques used to create 
the synthetic spectra of our simulations and discuss differ¬ 
ences between the different methods of estimating galaxy 
properties. In a companion paper (Guidi et ah, in prep., 
Paperll hereafter), we focus on the comparison of our simu¬ 
lated galaxies with data from the SDSS dataset. This paper 
is structured as follows. In Section 2 we describe the simu¬ 
lations and the simulated galaxy sample. Section 3 explains 
the techniques to create their SEDs and discusses uncertain¬ 
ties and biases affecting them, and in Section 4 we compare 
the magnitudes, stellar masses, gas/stellar metallicities, stel¬ 
lar ages and SFRs obtained using various methods. Finally, 
in Section 5, we give our conclusions. 


2 THE SIMULATIONS 

In this work, we use three sets of galaxy simulations con¬ 
sisting in total of fifteen galaxies formed in a ACDM uni¬ 
verse. Each set comprises the same five galaxies but adopts 
a different modelling of chemical enrichment and feedback, 
which is known to introduce differences in the final prop¬ 
erties of galaxies and in their evolution. The focus of this 
work is to identify whether the analysis techniques used to 
extract galaxy properties from simulations and observations 
introduces important biases that make the comparisons un¬ 
reliable. Our set of fifteen galaxies is ideal for this purpose 
as they all have similar total mass but span a wide range 
in gas/stellar metallicities, stellar ages and star formation 
rates. 

As the goal of this work is not to decide how realistic the 
galaxies are (which we discuss in Paperll), covering a variety 
of galaxy properties allows one to test the biases properly, 
unaffected by particular details of a given implementation. 
We use our sample to test minimum and maximum biases 
introduced in the conversion of simulations into observables, 
as this conversion is expected to primarily depend on the age 
of the stellar populations, the stellar masses, the metallicities 
and the galaxy morphologies. The latter is important as, on 
one side, dust will influence differently observations of face- 
on and edge-on galaxies and, on the other hand, the presence 
of gradients in galaxy properties might affect the derivation 
of their global properties from the observationally-obtained 
fiber quantities, since fiber spectrographs such as those used 
for SDSS sample only the inner region of galaxies. 

The initial conditions correspond to five galaxies which 
are the hydrody namical counterpart s of (a subset of) the 
Aquarius halos dSprinael et al.ll2008ll . These are similar in 
mass to the Milky Way and formed in isolated environments 
(no neighbour exceeding half their mass within 1.4 Mpc at 
redshift 2 = 0), but h ave different merger and a ccretion his¬ 
tories, as discussed in IScannapieco et al.l (l2009ll . The galax¬ 
ies have virial masses between 0.7 and 1.7 x IO^^Mq (cal¬ 
culated within the radius where the density contrast is 200 
times the critical density), stellar masses of 1 —10 x IO^^Mq, 
and gas masses of 3 — 10 x lO^^M©. 

For our first set of simulations, we have use d the ex¬ 
tende d version of the Tree-PM SPH code Gadget-3 (ISpringell 
l2005lf ■ which includes star formation, chemical enrichment, 
supernova Type la and Typell feedback, metal-dependent 
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cooling and a multiphase model for the gas component 
which allows the coexist e nce o f dense and diffuse phases 
l|ScannaDieco et ^ l200Sl. l200fil'l. These simu lations have 
been first p resented inIScan napiec o et alJ ll2 009l) a nd further 
analys ed in IScannapieco et al. ( 2010l ) and lScannapieco et alJ 
|20 i 3), and will be called throughout this paper A(-E)-CS 
or CS galaxies. The Scannapieco et al. code has been ex¬ 
tensively used for simulating galaxies of a wide range of to¬ 
tal/stellar masses, and shown to be successful in reproduc¬ 
ing the formation of galaxy discs from cosmological initial 
cond itions of Milky Way-mass ga l axies JScannapieco et al.l 
20081. see a l so ISawala et ^ l201( I: IScanna pieco et ah 20121;' 
Nuza et ^ l20l4 Creasev et al. l2015l : IScannapieco et al.l 
2015h . 


For our second simulation set, we used an updated 
version of the model of Scannapieco et ah, in relation to 
the treatment of chemical enrichment (Poulhazan et ah, 
in prep). These simulations will be referred to as A(-E)- 
CS^ or CS^. The updated code implements a different Ini¬ 
tial Mass Fun ction (Chabrier instead of Salpete r), chemi¬ 
cal yields from lPortmnnj_Chipsi_&_^ressM 1 19981 (while the 
CS model uses the lWooslev fc WeavenI 19951 yields), and the 
treatment of feedback from stars in the AGB phase, which 
contribute significant amounts of given chemical elements, 
such as carbon and nitrogen. For these reasons, the CS’'" 
galaxies have systematically higher chemical abundances 
compared to those in the CS sample. The modelling of en¬ 
ergy feedback is the same as in the standard Scannapieco et 
al. code. 


The final set of simulations, which will be referred 
to as A(-E) - MA or the MA galaxies, have used the 
lAumer e t al. (l2 013l'l code , whic h is an independent update 
to the Scannapieco et all (l2006l l model. This code has a dif¬ 
ferent set of chemical choices in relation to the initial mass 
function and chemical yields, includes also stars in the AGB 
phase, and assumes a different cooling function. More impor¬ 
tantly, it has a different treatment of supernova energy feed¬ 
back: unlike in the Scannapieco et al. code, where feedback 
is purely thermal, in Aumer et al. the supernova energy is 
diyided into a thermal and a kinetic part, and also includes 
feedback from radiation pressure. This results in stronger 
feedback effects compared to the Scannapieco et al. model 
and, as a consequence, produces galaxies that are more disk 
dominated, younger and more metal rich compared to the 
rest of our s i mulat ions. We refer the interested reader to 
lAumer et al. liiH) for full details on this implementation. 


The assumed cosmological parameters of the simula¬ 
tions are as follows: fim = 0.25, IIa = 0.75, fib = 0.04, 
as = 0.9 and Ho = 100/i km s“^Mpc“^, with h = 0.73. All 
simulations have similar mass resolution (2 — 5 x 10®Mq for 
stellar/gas particles and 1 — 2 x IO^Mq for dark matter par¬ 
ticles) and adopt similar gravitational softenings (300 — 700 
pc). 

Fig. [T] shows the color-composite images of the fifteen 
simulated galaxies, in face-on and edge-on views, obtained 
with the code SUNRISE (see below). The composite images 
in the (u, r, 2 :)-bands are genera ted using the algorithm de¬ 
scribed in iLupton et al.l J2004l 'l. The edge-on and face-on 
views are defined such that the total angular momentum of 
the stars in the galaxy is aligned with the z-direction. 


3 CREATING SYNTHETIC SPECTRA OF 

THE SIMULATED GALAXIES 

To derive the observables of our simulated galaxies and in¬ 
vestigate the biases introduced in the process, we follow 
three approaches. First, we use different Stellar Population 
Synthesis (SPS) models, which provide the resulting spec¬ 
trum due to the emission of the stars, as well as differ¬ 
ent informations about the stellar populations at a given 
metallicity and age (e.g. mass loss, ionizing photon flux, 
Lick indices, number of black holes and neutron stars). Sec¬ 
ond, we add a simple dust model to the prediction of the 
SPSs, in order to better compare the resulting spectra from 
those of observed galaxies. Finally, we postprocess the sim¬ 
ulations with the radiatiye transfer code sunrise djonssonl 
I 2 OO 6 I : [jonsson. Groves fc Coxll2009l l. which is computation¬ 
ally slower but more consistent with the underlying hydrody- 
namical simulation. This gives the full SED including stellar 
absorption features, nebular emission, and extinction due to 
dust as light travels out through the interstellar medium. 

The SPS models are commonly used as a postpro¬ 
cessing of hydrodynamical simulations as they provides 
a fast and easy way to estimate the light distribution 
of an ensemble of stars. The use of radiative transfer 
codes has also been applied to hydrodynamical simula¬ 
tions (e.g. iGover nato et al]|2009l: IScanna pieco et all 120101 : 


iHayward et al. 20131 : Christensen et ahT 20141 '). This ap¬ 

proach has the advantage that can consistently consider the 
distribution of dust (traced by the metals) although it also 
needs to introduce certain assumptions and simplifications. 

In the following subsections, we describe the main as¬ 
pects of these models and the assumptions that are relevant 
for our work, as well as discussing the main uncertainties 
and sources of biases in the creation of the synthetic SEDs. 


3.1 Stellar population synthesis Models 

We use various SPS models to create synthetic spectra of 
our simulated galaxies. The spectra are obtained assuming 
that each star particle represents a Simple Stellar Popula¬ 
tion (SSP) parametrized by the age and metallicity of the 
star, and normalized by the mass of the particle. In the 
simulations, the mass of a star particle changes with time 
owing to the supernova ejecta and the mass loss during the 
AGB phase; the mass normalization should be made with 
the mass of the stars at the formation time (note that this 
is important in order to avoid double-counting the effects of 
stellar mass loss). The spectrum of the galaxy is obtained 
by summing up the spectra of individual stars. The galaxy 
spectra can be convolved with given photometric bands and 
integrated to get the magnitude of the galaxy in the bands. 

We use 5 different SPS models to obtain spectra, mag¬ 
nitudes and colors of our simulated galaxies, which allows us 
to identify systematic effects due to their varying assump¬ 
tions. Although moderate, some differences between these 
models are expected, as uncertainties in modelling the stel¬ 
lar evolution still exist (e.g. different treatment of convec¬ 
tion, rotation, mass-loss, thermal pulses during AGB evo¬ 
lution, close binary interaction), as well as in the empirical 
and theoretical stel l ar spectral lib raries (see the reviews by 
IWalcher et al.l[201ll : IConrovll2013l l. The main characteristics 
of the SPS models are described in the following, and in Ta- 
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A-CS A-CS+ A-MA 



D-CS D-CS+ D-MA 



E-CS E-CS+ E-MA 



Figure 1. Multi-band {u,r,z) images of our fifteen galaxies, for face-on and edge-on views. 


ble 1 we give a summary of the input parameters we chose, 
taken as homogeneous as possible to make the interpretation 
of our results more clear. 

• BC03 llBruzual fc Charlotj l200lj ') computes, using dif¬ 

ferent evolutionary tracks, the spectral evolution of stellar 
populations at a resolution of 3 4 FWHM in the optical and 
at lower resolution in the full wavel ength range. _ 

• SB99 (“STARBURS T99”, iLeitherer et al.l 1 19991: 

I Vazquez fc LeithereJ l2005l l is a web-based platform that 
allows users to run customized SPS models for a wide 
range of IMF, isochrones, model atmospheres, ages and 
metallicities. The spectral resolution reaches ~ 1 A FWHM 
in the optical wavelength range for the fully theoretical 
spectra. _ 

• PE f“PEGASE”. rFioc fc Rocca-Volmerangdll997l . ll999l l 
uses an algorithm that can accurately follow the stellar 


tracks of very rapid evolutionary phases such as red super¬ 
giants or TP-AGB. The stellar library includes also cold star 
parameters, and a simple model for nebular emission (con¬ 
tinuum + lines) can be added to the stellar spectrum. 

• FSPS (“Flexible S tellar Population Synthesis 2.3”, 
IConrov. Gunn fc Whitell2009l l is a flexible SPS package that 
can compute spectra at resolving power A/AA « 200 — 500. 
In addition to the choice of IMF, metallicities, ages, the user 
can select a variety of assumptions on Horizontal Branch 
(HB) morphology, blue straggler population, location of the 
TP-AGB phase in the HR-diagram and post-AGB phase. 

• M05 llMarastonI l2005l l is different from the other 
SPS models in the treatment of Post Main Sequence 
stars, namely using t he Fuel Consumption theorem 
llRenzini fc BuzzonilllOSd l to evaluate the energetics. It gives 
spectra at a resolution 5 — 10 A in the visual region and at 
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20 — 100 A from the NUV to the near-IR for either blue or 
red populations on the horizontal branch. 


3.2 Dust 

Dust extinction is an important ingredient in the estimation 
of observables from the simulations, as in observed galaxies 
dust effects can be large, especially for edge-on systems. Dif¬ 
ferent authors ha ve modelled d ust extinction curves of the 
Milky Way (e.g. ISeatonI Il979l : ICardelli. Clavton fc Matlii^ 
1989(1. of the Large and Small Ma gellanic clouds (e.g. 


FitzpatrickI 1 19861 : [Bouche^j^^ 1985ll. and also of exter¬ 


nal galaxies (e.g. Silva et al.l [19981: IC harlot fc Falll I 2 OOOI : 


ICalzetti et al.l 200CI : Fischera fc Dopita 2005l l. which are of¬ 
ten used to correct for dust extinction. 

In our calculations of the dust -corrected mag n itudes 
and colors, we use the model of ICharlot fc Falll ll2000l . 
C FOO hereafter) in the slightly differ ent formulation given 
in Ida Cunha. Chariot fc Elbaa (l2008l . dC08 hereafter), that 
consistently extend the CFOO model to include dust emis¬ 
sion, in order to allow the interpretation of the full UV-far 
infrared galaxy SEDs. 

CFOO is an angle-averaged time-dependent model, with 
extinction curve depending on the wavelength, and on 
the Stellar Population (SP) age. This model has been 
proven to work reasonably well for a wide class of galax- 
ies, and it has been alre ady implemented in SPS models 
dBruzual fc Charl^ l2003l 'l. In CFOO, stars are assumed to 
born in “birth c louds” that disp erse after a certain amount 
of time (see also lSilva et al.lfl998h : the transmission function 
of the SP is the product of the transmission function in the 
birth cloud (which depends on the SP age t) and in the ISM: 

^ I exp [-(1 - n) TV ( 55 ^)“'^'*"’] for t ^ tec 
(1 for t > Ibc 


rrnlSM 

T\ = exp 


-^.Tv 


A 


5500 A 


—Nism 


where Ibc is the birth cloud life-time (see CFOO and dC08 
for details). The resulting attenuated luminosity Lx given 
the intrinsic luminosity S\ is then 

Lx{t) = Sa 

We set the free parameters to the values given in dC08 for 
normal star-forming galaxies, namely: tnc = 10 Myr, tv = 
1.5, = 0.3, A^bc = 1-3 and Nism = 0.7. 


3.3 Radiative Transfer 

To calculate the full far-UV to submillimeter SED of our 
simulated galaxies we use the Monte Carlo Radiative Trans- 
fer code sunrise in the post-processing p hase. SUNRISE 
(|jonssonll2006l : [j onsson. Groves fc Coxll2009h is a 3D adap¬ 
tive grid polychromatic Monte Carlo radiative transfer code, 
suited to process hydrodynamical simulations, sunrise as¬ 
signs a spectrum to each star particle in the simulation, 
and then propagates photon “packets” from these sources 
through the dusty ISM using a Monte Carlo approach, as¬ 
suming a constant dust-to - meta ls mass ratio, which we fixed 
to 0.4 according to iDwe 5 (1 19981') . 


In the standard sunrise implementation, each stellar 
particle older than 10 Myr is assigned a spectrum corre¬ 
sponding to its age and metallicity from the input SPS 
model, in our case SB99 (see Table 1). Star particles younger 
than 10 Myr are assumed to be located in their birth 
clouds of molecular gas, and are given a modified spec¬ 
trum which accounts for the effects of HII and photo¬ 
dissociation regions (PDRs). The evolution of HII regions 
and PDRs are described by the photo-ionization code map¬ 
ping s III llCroves. Dopita fc Sutherlandl[^04l : ICroves et all 
l2008li . which is used to calculate the propagation of the 
source spectrum through its nebula. The HII regions ab¬ 
sorb effectively all ionizing radiation and are the sources of 
hydrogen recombination lines, as well as hot-dust emission. 
The only mappings III parameter not constrained by the 
hydrodynamical simulation is the PDR clearing timescale; 
in sunrise this free parameter has been changed to the time- 
averaged fraction of stellar cluster solid angle covered by 
the PDR (/pdr) , for which we use the fiduci dal value of 
/pDR = 0.2 as in ijonsson. Groves fc CoxI (l2009l ~). 

Dust extinction in sunrise is described by a Milky 
Way-like extinction curve normalized to R y = 3.1 

dCardelli. Clavton fc Mathij 1 19891 : IPraind I 2 OO 3 I I which in- 
clude s also the 2175A bump observed in our galaxy lICzernvl 
1200711 . Model cameras are placed around the simulated 
galaxies, to sample a range of viewing angles. The emer¬ 
gent flux is determined by the number of photons that flow 
from the galaxy unscattered in a given camera’s direction, 
as well as those scattered into the line of sight or reemitted 
in the infrared by dust into the camera. In the calculations 
presented here we use two cameras, one looking face-on and 
the other edge-on to the galaxy, to take into account the two 
extreme cases in terms of optical depth. Finally, the flux in 
the cameras is convolved with the bandpass filters and inte¬ 
grated to get the broadband images and magnitudes in the 
chosen photometric bands. 

Fig. [2] is an example of the spectra obtained for one of 
our simulated galaxies, C-MA, showing in the upper panel 
the SB99 spectrum both without dust (pink line) and cor¬ 
rected with the CFOO model (blue line), as well as the dust- 
free stellar spectrum with sunrise (black dashed line). The 
lower panel shows again the sunrise stellar spectrum with¬ 
out dust, the full nebular-l-stellar spectrum without dust 
(red), and the full spectrum with dust from the edge-on 
camera (green). 


3.4 Uncertainties and biases 

In addition to the different assumptions of the SPS mod¬ 
els and the treatment of dust and radiative transfer, there 
are a number of choices on the process of constructing the 
synthetic SEDs of simulated galaxies. In some cases, these 
constitute important sources of biases that need to be kept 
in mind when results are interpreted, while others do not 
strongly affect the final SEDs. When appropriate we will 
test and quantify these effects, which we describe below. 


3.4-1 Aperture bias 

While simulations give the full information on the phase 
space of particles that constitute a galaxy, observations 



































































6 Guidi, Scannapieco & Walcher 


Table 1. Summary of the characteristics and choices of parameters for the 5 SPS models used in our work. The input stellar model in 
SUNRISE is SB99 with the same parameters. 


Model 

IMF 

Age range (yr) 

Metallicity range 

Wavelength range 

Stellar tracks 

Stellar library 

BC03 

M05 

SB99 

PE 

FSPS 

Chabrier^^^ 

Kroupa^^^ 

Kroupa^^^ 

Kroupa^^^ 

Kroupa^^^ 

10^ - 2 ■ 10^“ 

10^ - 1.5 ■ 
lO'* - 1.5 ■ IQio 

0 - 2 ■ 10i° 

3 ■ lO'^ - 1.5 ■ 10^“ 

0.0001 - 0.05 

0.0001 - 0.07 

0.0004 - 0.05 

0.0001 - 0.05 

0.0002 - 0.03 

91 — 160 fim 

91 — 160 fim 

91 — 160 fim 

91 A — 160 fim 

91 A — 10 mm 

Padova 1994 
Cassisi / Schaller 
Padova 1994 

Padova 1994 

Padova 2007 

BaSeLS.l / STELIB 

BaSeL3.1 / Langon & Mouchine 
Pauldrach / Hillier 

BaSeLS.l / ELODIE 

BaSeLS.l / Langon Sz Mouchine 

NOTES: 

(1) Mass range: 

m = 0.1-100 M© , 

a = 2.3 for m > 1 Mq; (2) a = 1.3 for 

m =0.1-0.5 M0, 

a = 2.3 for m = 0.5-100 Mq. 



3000 4000 5000 6000 7000 8000 9000 10000 

Wavelength [A] 


Figure 2. Examples of the SEDs obtained for one of our simulated galaxies (C-MA) using different methods, including the stellar 
spectrum only (SB99), the stellar spectrum and the CFOO simple dust model (SB99 + CFOO), the stellar SUNRISE dust-free spectrum 
(SRnoiSM stellar) and total SUNRISE spectra in the absence (SRhoISm) uud presence of dust, for the edge-on view (SRedgeon)- The lower 
panel shows the SDSS filter transmission functions for the u,g,r,i and 2 -bands. 


usually have limitations related to the region that can be 
observed/measured, and the different techniques used to 
define this region in obser vations and simulati ons can in¬ 
troduce some systematics (IStevens et al.l[2014l) . In SDSS, 
the properties derived from spectra (e.g. stellar age and 
metallicity, gas metallicity, SFRs) are affected by the 
small aperture (3” arcsec) of the SDSS spectrograph, 
which samples only the inner part of the target galax¬ 
ies. The presen ce of metallicity gradients in galaxy proper- 
ties and colors jlBeU^^eJonjlgOOOliJPiUd ^ et al.l[2012l : 


IWelikala fc Knei? 20121 :' Sanchez et al. 2015 ) can therefore 

lead to substantial uncertainties in the observed measure¬ 
ments when one wants to estimate the total quantities. The 
implic ations of this bias have been discussed by several au¬ 


thors (iBell fc de Jong||2000l: KochanekjPahre_^_^hlc^[200^ 


Bajdry_et_ajjj200j|_ [C^mezet al.l l2003l : Brinchmann et ahl 
20041 : MacArthur et al. 20041 . e.g.); in some cases methods 
for aperture correction which exploit the spatially-resolved 
color informatio n can be used to extrapo l ate to global quan - 
tities (e.g. SFRs. lBrinchmann et al.l2004l . [^alim et al.l200lt) . 


We used our simulations to test the effects of using a 
single fiber in the derivation of properties such as metallic- 
ities and stellar ages, by including in the calculations both 
all particles in a 60x60 kpc field of view (FoV) in the face- 


on projection (full FoV), or only those within a small region 
that mimics the size sampled by the SDSS fiber spectro¬ 
graph, i.e. a circular region of 4 kpc radius with the galaxy 
in the center observed face-on (fiber FoV (Qi which according 
to the cosmology adopted in our simulations corresponds to 
redshift 2 ~ 0.15 for the SDSS instrument aperture. Note 
that we are limited to further decrease the size of the fiber 
field of view by having enough particles to reach a good 
statistical sample, in particular in the case of gas particles. 

Figs. El m and El show the number-density histograms 
of stellar ages and stellar/gas metallicities for the 15 simu¬ 
lated galaxies, when we consider all particles in the galaxy 
and only those within the fiber. For each galaxy, we also 
show the fraction between the number of star/gas particles 
within the fiber to the total number of star/gas particles in 
each bin, nfib/ji-tot, that we refer to as the sampling func¬ 
tion. In the case of stellar ages we find that, in general, 
fiber quantities preferentially sample the older populations, 
with sampling functions typically higher than 50%. In con¬ 
trast, the young populations are sampled in various ways, 
depending on the age profile of the galaxy. In some cases. 


^ In the z-direction we count all stars in the halo as identified 
with SUNFIND as bound to the halo. 
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Figure 3. Number-density histograms of stellar ages for the 15 simulated galaxies, considering all star particles in the galaxies as well as 
only those within the fiber. For each galaxy, we also show the fraction between the number of stars within the fiber to the total number 
of stars in each bin, ngb/ntot, that we refer to as the sampling function. 


young populations are not at all sampled and will be com¬ 
pletely missed in the calculations using particles within the 
fiber. It is clear, in any case, that the sampling within the 
fiber is not constant or even similar for all galaxies, which 
means that it is not possible to reliable estimate the mean 
stellar age of a galaxy from the fiber quantities, or even es¬ 
timate the fiber bias without having additional information 
(e.g. importance of age/metallicity profiles) of the regions 
outside the fiber. We note that the relative contribution of 
old and young populations in a galaxy will in general reflect 
the relative importance of bulges and disks, but with fiber 
quantities this information might not be a reliable reflec¬ 
tion of the real relative contribution of the different stellar 
components of a galaxy. 


pared to those in the MA sample, with the CS”'" galaxies 
lying in between. This results from the different implemen¬ 
tation of chemical enrichment (chemical yields increase from 
the CS galaxies, to CS^ and then to MA), chemical diffu¬ 
sion (only included in the MA simulations) and feedback 
(from weaker in CS to stronger in MA, with CS^ in be¬ 
tween). In particular, the MA galaxies have more peaked 
dist ributions, and have less strong metallicity gradients (see 
also lAum er et al.||2013l) compared to the CS and CS”*" galax¬ 
ies (see also Tissera. White fc Scannapiecdl2012l l. Unlike for 
the stellar ages, for the stellar metallicities we find that the 
sampling shows much less variation from galaxy to galaxy, 
as in all cases there is a preferential sampling of the more 
metal-rich populations. 


Similar considerations can be made in the case of the 
stellar metallicities, shown in Fig. [d] The metallicity distri¬ 
butions of the galaxies in the CS sample is broader com¬ 


The distributions of gas metallicities are more complex, 
and show important differences depending on the details of 
the implementation of chemical enrichment (IMF, chemical 
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Figure 4. Number-density histograms of stellar metallicity (in solar units, assuming Zq = 0.02), both for all star particles and inside 
the fiber field of view. For each galaxy, we also show the fraction between the number of stars within the fiber to the total number of 
stars in each bin, riflu/utot, that we refer to as the sampling function. 


yields, etc.) and on the absence/presence of chemical diffn- 
sion, which is included only in the MA sample and leads 
to a much less broad distribution of gas metallicities (note 
the different x-scales in the lower panels). From Fig. [5l it 
is clear that the sampling of gas metallicities varies signif¬ 
icantly from galaxy to galaxy. While in general there is a 
preferential sampling of the metal-rich regions, gas parti¬ 
cles with intermediate metallicities can also contribute sig¬ 
nificantly to the final metallicity when fiber quantities are 
considered. For the MA sample, we find that, within the 
fiber, only a narrow range of gas metallicities are sampled, 
which reflects the effects of metal diffusion that tends to 
give a smoother metallicity distribution. In the case of the 
D-CS galaxy, we also find that although the metallicity dis¬ 
tribution is broad, only a very narrow range of metallicities 
are sampled within the fiber. According to our findings, for 
estimating the gas metallicy of a galaxy it is of primary im¬ 


portance that the fiber bias can be quantified; otherwise it is 
not possible to derive the mean gas metallicity of the whole 
galaxy reliably. 

In summary, our results show that the bias due to 
the fiber is in general reflected in the tendency to sample 
older and more metal-rich stellar populations and metal- 
enriched regions of the ISM, although the shape of the sam¬ 
pling function highly depends on the significance of gradients 
in ages/metallicities (which also somewhat depend on the 
modelling of physical processes in the simulations), affect¬ 
ing in particular galaxies with stronger gradients. 


3 . 4.2 Variation of the assumed Initial Mass Function 

The choice of the Initial Mass Function (IMF) can 

strongly bias the de r ivation of many observab les _ of 

galaxies (see ISalpeterl Il955l : iMiller fc Scald 1 19791 : IScald 
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Figure 5. Number-density histograms of gas metallicity, both for all star particles and inside the fiber field of view. For each galaxy, we 
also show the fraction between the number of stars within the fiber to the total number of stars in each bin, that we refer to 

as the sampling function. Note the different x-range plotted for the CS/CS"^ and MA samples. 


1 19861: H roupal l2002l ^. for example, the assumed IMF 
has a direct influence on the M/L ratio depending on 
the age, composition and past star for mation history 
llBaldrv fc Glazebrookll200^ : IChabrierll2003l 'l. which in turn 
influences the stellar mass derivation. The IMF is mea¬ 
sured in many environments, still there is not yet consen¬ 
sus on a number of impor tant aspects, such as whether 
it is universal or not (e.g. Hoversten fc GlazebrookI l2008l : 


iBastian, Covey fc Meveij|2oid'l.' 


We assume in general a Kroupa IMF with slope a = 1.3 
for masses 0.1-0.5 Mq and a = 2.3 in the mass range 0.5- 
100 Mq when we apply the different SPS models, except 
in the case of BC03, where we instead use a Chabrier IMF 
(Table 1). The difference between these two IMFs is mod¬ 
erate, however, due to their very sim ilar shape, particularly 
at the high-mass end (see Fig. 17 in lEllisl[2008l 'l. Note that 
our three galaxy samples use different IMFs to calculate 


feedback and chemical enrichment: Salpeter (CS), Chabrier 
(CS”'") and Kroupa (MA). 


4 RESULTS 

In this section we describe the methods we use to derive the 
galaxy properties - magnitudes and colors, stellar masses, 
stellar ages, stellar and gas metallicities and star forma¬ 
tion rates - from their synthetic spectra, and make a de¬ 
tailed comparison between results obtained with the differ¬ 
ent methods, in order to identify biases and systematics. For 
each property, different techniques are used in the deriva¬ 
tion, the selection of which has been made taking into ac¬ 
count the questions we want to answer: (i) what is the range 
of variation in observationally-derived quantities when dif¬ 
ferent techniques are used, (ii) do they agree with the ’real’ 
quantity; and (Hi) how do observational definitions of galaxy 
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properties effect their measured values. In the case of obser¬ 
vations, we particularly focus on the techniques employed in 
SDSS, which will be used in Paperll to compare simulated 
and observed galaxies. 

In order to avoid additional biases, we always use the 
same field of view defined as a 60 kpcxGO kpc region with 
the galaxy in the center both when we derive the properties 
directly from the simulations and from the SPS/sunrise 
spectra. 

4.1 Magnitudes and colors 

We first compare the colors and the absolute magnitudes of 
our simulated galaxies, in the 5 SDSS photometric bands 
(u,g,r,i,z, see Fig. obtained using different methods, 
namely: 

• BC03, SB99, FSPS, PE, M05: these refer to the 
magnitudes calculated by applying the five different SPS 
models described in Section lSTTI As explained in the previous 
section, for each star particle the magnitudes are obtained 
via a linear interpolation of the SPS tables according to its 
age and metallicity, normalizing with the particle mass at 
the formation time. 

• CFOO: we include dust effects to BC03 using the model 
of CF00/dC08 described in Section [3.21 

• SRnoisM, SRfaceon) SRedgeon: results from the radiative 
transfer code SUNRISE as explained in Section 13.31 in the 
absence of dusl[3, and for the face-on and edge-on views in¬ 
cluding dust, respectively. Note that the input SPS model 
for SUNRISE is SB99. 

The left-hand panel of Fig. compares the magnitudes 
obtained using the different SPS methods and sunrise, 
when we ignore the effects of dust. We show the differences 
with respect to BCOi0 as a function of the stellar metallic¬ 
ity of the galaxj0 (in solar units, assuming Zq = 0.02), as 
systematics in SPS models are expected to increase both at 
low/high metallicity (as well as at younger stellar age s, e.g., 

I Conroy, Gunn fc White! l2009l : I Conroy fc GunnI l2010l ~) . The 
different SPS models show in general very good agreement, 
with magnitude differences of 0.1. We detect however 
some systematics: models FSPS and M05 usually predict 
lower magnitudes (i.e. brighter galaxies) compared to BC03, 
while models SB99 and PE give systematically higher mag¬ 
nitudes than BC03. As expected, differences are somewhat 
larger (but still moderate) for the most metal-rich galax¬ 
ies, which are also those that exhibit a higher fraction of 
young (ages < 10 Myr) and intermediate (ages in the range 
[0.1 — 2] Gyr) stellar populations. For these ages the un¬ 
certainties in the treatment of young stars, post red-giant 
branch and TP-AGB phase stars are larger. The galaxy 

^ We note, however, that dust is included through the sub-grid 
dust model around young stars in MAPPINGS. 

® We note that, even if BC03 assumes a different IMF (Chabrier) 
compared to the other SPS models, this does not introduce any 
significant systematic effect in the derived magnitudes/colors, 
since for the more luminous stars (M* > IMq) the IMF slope 
is the same as Kroupa (see Table 1). 

^ This is calculated as the average mass-weighted metallicity over 
all stellar particles in the simulated galaxy (method SIM in Sec¬ 
tion 14.311 . 


which exhibits the largest AM (in all bands) is A-MA, which 
has log{Z/ZQ) « —0.26 the most extreme SFR f'Section ld.SI) 
and the youngest mean stellar age in our sample (see next 
sections). 

The results of SRnoiSM also agree well with the rest of 
the models, particularly with those of SB99, as SB99 is the 
input SPS model of sunrise. The remaining differences are 
because in sunrise the spectrum for young stars is calcu¬ 
lated with MAPPINGS III, which includes the contribution 
from nebular continuum and emission as well as dust absorp¬ 
tion and IR emission that is modeled as sub-grid physics. 
The largest differences are found for A-MA, the galaxy with 
the highest number of young star particles. In fact, on one 
side the ionizing photons are reprocessed as nebular emission 
lines and continuum increasing the total flux in the optical 
and, on the other hand, the young particles are also more 
extincted due to the sub-grid treatment of dust absorption 
in MAPPINGS III, giving total magnitudes in general lower 
(brighter) by < 0.1 mag than modelling the SED including 
only stellar emission. 

The right-hand panel of Fig. [5] shows our results in 
the case of including the effects of dust, either with the 
CF00/dC08 correction or, more consistently, with the SUN¬ 
RISE code. For comparison, the magnitude differences are 
calculated with respect to the BC03 model. When the CFOO 
simple dust model is included we find, as expected, fainter 
galaxies, with differences between ~ 0.6 — 0.8 mag for the 
u and g-bands and of ~ 0.4 — 0.5 mag for i, r and 2 . As 
the CFOO model is angle-averaged, there is no difference be¬ 
tween face-on and edge-on views. When we consider radia¬ 
tive transfer effects in SUNRISE, we detect almost no differ¬ 
ence if we see the galaxies face-on (SRfaceon) or if dust is 
ignored (SRnoiSM method, included here for comparison). 

As expected, larger differences are detected when galax¬ 
ies are seen edge-on, as dust effects are maximal in this 
case. However, we find significant differences only for 
metal-rich galaxies, while for galaxies that are metal-poor 
(log(Z/ZQ) < -0.45) dust effects are unimportant. This is 
because the amount of dust in SUNRISE is directly propor¬ 
tional to the gas metallicity (galaxies with low stellar metal¬ 
licity also have low gas metallicity). Observations also con¬ 
firm that the dust-to-metals ratio is nearly constant over 
a large range of metallicities and r edshifts, and dust ef¬ 
fects are small in met al-poor galaxies dZafar fc Watsonl2013l . 
iMattsson et al.|[2014l') . For the metal-rich galaxies, the inclu¬ 
sion of dust makes them fainter by factors of ~ 0.8 — 1 dex 
in the u-band and of 0.3 — 0.5 dex for the 2 -band. In¬ 
terestingly, for the most metal-rich galaxies the magnitudes 
derived from the SRedgeon and CFOO methods are very sim¬ 
ilar. Note also that the CS”*" galaxies have a systematically 
lower reddenning compared to the MA ones, even at high 
metallicity. This is explained by the lower amount of met¬ 
als in the ISM, and therefore lower amount of dust of the 
galaxies in the CS”'" sample compared to those in MA. 

To explore the dependence of dust effects on metallicity 
for the metal-poor sample, we rerun sunrise for one galaxy 
(C-CS), using a higher dust-to-metals ratio of 4 (10 times 
larger than t he standard value of 0.4 assumed in our sunrise 
calculations. lDwe^ll998l l. The results are shown as open and 
filled squares, respectively for the edge-on and face-on views. 
In this case, we indeed find little effects in the case of the 
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Figure 6. Difference between the absolute magnitudes of the simulated galaxies obtained using different methods and those predicted 
by the dust-free BC03 model, in the 5 SDSS bands and as a function of stellar metallicity. The left-hand panel shows results for dust-free 
models and the right-hand panel compares results of the CFOO simple dust model and of the SUNRISE code. The open/filled black squares 
are results for C-CS, when a higher dust-to-metals ratio is assumed in the SUNRISE calculations in edge-on and face-on views, respectively. 
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Figure 7. Difference between the colors of the simulated galaxies obtained using different methods and those predicted by the dust-free 
BC03 model, as a function of stellar metallicity. The left-hand panel shows results for dust-free models and the right-hand panel compares 
results of the CFOO dust model and of the SUNRISE code. 
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face-on view and larger differences in the edge-on case, with 
results very similar to those of CFOO. 

The differences in magnitndes, when different methods 
are applied, are translated into differences in the colors, as 
shown in Fig. [3 Both in the absence (left-hand panel) and 
in the presence (right-hand panel) of dust, differences are 
larger for the blner colors {u — g, u — r and g — r), and for 
the most metal-rich galaxies. 

Our results show that applying different SPS to sim¬ 
ulations gives visual magnitudes of galaxies with spread 
0.05 — 0.1 dex for metal-poor galaxies and of 0.1 — 0.2 for 
metal-rich ones, depending on the model, with BC03 ap¬ 
pearing as intermediate among the 5 SPS models tested. 
Furthermore, while for face-on galaxies or edge-on galaxies 
with low metallicities {\og{Z/Z q) < — 0.45) the effects of 
dust can he in first approximation ignored, when galaxies are 
seen edge-on and have a significant amount of dust (i.e. are 
metal-rich), errors in the magnitudes can be up to 0.7 — 1.5 
dex if dust is ignored, depending on the band. This means 
that for edge-on galaxies, particularly metal-rich ones, re¬ 
liable magnitudes will not be possible to obtain without a 
proper modelling of the dust extinction. Similar considera¬ 
tions can be made for the colors, that can not be calculated 
with precision better than 0.1 — 0.5 dex. 

The previous figures showed the dependence of the 
galaxies’ magnitudes on different models and assumptions. 
In simulation studies, it is common to use SPS models to 
convert masses into luminosities, and to compare models 
with observational data (e.g. light proHles, magnitudes, etc). 
We have shown that the use of different SPS models intro¬ 
duces some changes in the predicted magnitudes, that are 
however moderate. On the other hand, when dust effects are 
included or a proper radiative transfer treatment is consid¬ 
ered, larger differences might appear. In the left-hand panel 
of Fig. [8] we show the position of our simulated galaxies 
in the color-magnitude diagram; in this case in the Mr vs 
[u — r) plane. We show the results for BC03, the spread 
found for all SPS models in both plotted quantities, and the 
results for the BC03 model with the CFOO correction. We 
also show results for SRnoiSM, SRfaceon and SRedgeon- These 
cover the commonly used ways to post-process simulation 
results, and allow to understand what are the possible off¬ 
sets expected when more sophisticated methods are used to 
calculate the magnitudes. Most of our simulated galaxies 
have Mr between —22 and —19, and (u — r) color between 
1.5 (except for A-MA) and 2.3 (see also Paperll for a 
comparison with SDSS data). The A-MA galaxy lies at a dif¬ 
ferent position compared to the other galaxies, particularly 
when dust effects are ignored or when it is included but the 
galaxy is seen face-on. Note that this is a young galaxy with 
strong emission lines, and the Ha line falls in the r-filter at 
z = 0. A-MA however moves photometrically closer to the 
other galaxies if dust extinction is included. 

In observations, measuring the total magnitudes of 
galaxies can be plagued by different observational problems 
(e.g. low signal-to-noise, sky brightness, bad calibrations and 
sky subtraction), which affect particularly the outer part of 
the profiles where the S/N is lower. For this reason, def¬ 
initions on how to measure the total light are used in all 
galaxy surveys, although they ca n differ substantially among 
auth o rs and collaborations (e.g. iKronI Il98d . iBlanton et al.l 
I 2001 I . iNorberg et al.l[2002l l. We have calculated the mag¬ 


nitudes and colors of our simulated galaxies following the 
techniques used in SDSS, in particular ca lculating the so- 
called Petrosian an d Model magnitude s (IPetrosianI Il976l: 
Blanton_et_alJ_ 200l|^ Yasuda et al.l 1200ll : iKauffmann et aD 


2 OO 3 I : ^lim et alT 2007 1. These were derived from the im¬ 


ages obtained with sunrise, i ncluding dust and in the face- 
on view (SRfaceon), as follows l|Blanton et al.ll200ll 'l: 


• P ETRO: we calculate the Petrosian Radius (IPetrosiarJ 

in the r-band, and we derive the Petrosian magnitudes 
in all bands taking the flux inside Np = 2 Petrosian radii; 
as in SDSS, we assume a Petrosian Ratio Rp = 0.2. Note 
that, in general, this method samples only part of the total 
flux, the fraction d epending on the luminosity profile (see 
iGraham et ahllioOSl for detailed calculations). 

• MODEL: in this case, each image is matched to differ¬ 
ent luminosity profiles, and the magnitudes are calculated 
from the prohle which give s the best fit. We used the code 
GALFIT to perform the ht (IPeng et al.ll20ldl . assuming ar¬ 
bitrary axis ratio and position angle, and used two different 
profiles: Exponential and DeVaucouleur. Depending on the 
galaxy, one of these proHles provided the best ht: for 5 galax¬ 
ies this was an exponential, for 10 a DeVaucouleur prohle. 


The right-hand panel of Fig. |S]compares the results ob¬ 
tained with the PETRO and MODEL methods to those of 
SRfaceon. We hnd in general similar values of magnitudes 
for the PETRO and SRfaceon methods, with PETRO magni¬ 
tudes always higher (i.e. fainter) than those of SRfaceon as 
expected. In the case of the MODEL method, we also hnd 
fainter galaxies compared to the results of SRfaceon- In both 
cases, differences are lower than 0.3 dex for 9/15 galaxies, 
with a maximum of 0.6 dex for the remaining 6 systems. 
Colors are similar in the three models, with diherences in 
{u — r) < 0.1 — 0.2 dex. Our results show that, for our sim¬ 
ulated galaxies the PETRO and MODEL methods can not 
recover the real magnitudes/colors of galaxies with a preci¬ 
sion better than 0.2 — 0.3 dex. 


4.2 Stellar mass 

The stellar mass of galaxies is an important proxy of how the 
galax y populations evolved over cosmic time fe.g. lBell et al.l 
I 2 OO 3 II . In this section, we compare the stellar mass of our 
simulated galaxies, obtained in different ways, including the 
simple sum over the mass of stellar particles (as done in 
simulation studies) and those obtained using different post¬ 
processing techniques that mimic observations, in particular 
the ones used in SDSS. These are described below: 

• SIM: the total mass of star particles in the simulated 
galaxy. We include all stars in the 60 kpcx60 kpc held of 
view (see sec. 0, i.e., the same Held of view of sunrise. 

• BC03: the Hnal mass of star particles is calculated with 
the BC03 model, which considers the mass lost by a stellar 
population since it was formed, and until the present time. 
The BC03 Hnal mass of each star particle is obtained nor¬ 
malizing with its initial mass. The total stellar mass of the 
galaxy is the sum over the Hnal masses of the star particles 
within the Held of view. 

• SED: the total stellar mass is estimated by Htting the 
SUNRISE face-on optical galaxy spectrum (SRfaceon method) 
over the range 3800 — 9000 A using the starlight code 
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Figure 8. Color-magnitude diagram of the simulated galaxies, obtained using different methods to derive the synthetic spectra. The 
left-hand panel show results for the SPS, CFOO and SUNRISE methods and the right-hand panel compares results obtained following the 
SDSS techniques. 


dCid Fernandes et al.ll200^ ) with the BC03 SPS model. The 
spectra include nebular emission (masked during the fit) and 
are dust-extincted. 

• PHOTO: we fit the photometric (u, g, r, i, z)-band mag¬ 
nitudes obtained from the SRfaceon spectrum to a grid of 
models as updated from BC03 in 2007 (CB07, Chariot & 
Bruzual, priv. comm.); the grid spans a large range in 
galaxy star formati on histories, ages and metallicities (see 
IWalcher et al.lli^OOSl for a description of the code). To remove 
the nebular contribution from the broad-band magnitudes 
(the fitted model considers only stellar light), we calculate 
the relative contribution of nebular emission within the fiber 
in each photometric band fitting the fiber spectrum with the 
STARLIGHT code, and assume that the relative contribution 
of nebular emission for the total galaxy is the same as in the 
fiber. This procedure allows one to mimic the stellar mass 
estimation of the Garching SDSS DR7 (which we will use in 
Paperll to compare simulated and observed galaxies). 

• PETRO: Petrosian stellar masses are derived as in 
PHOTO-MASS, but using the Petrosian magnitudes (Sec- 
tion l4.1l) as input for the fit (i.e. within 2 times the Petrosian 
radius). The same procedure to remove nebular emission as 
in PHOTO-MASS was applied in this case. 

• MODEL: The same as PETRO-MASS, but using the 
Model magnitudes fSection 14.HI . 

In order to assess the effects of observational biases on 
the determination of the stellar mass of galaxies, it is impor¬ 
tant to consider the effects of mass loss of stellar particles 
in the simulations, as for all observationally-derived quan¬ 
tities, the fitting of spectra/magnitudes is done using the 
SPS models that include a mass loss prescription. Our 15 
simulations cover the cases where mass loss is not followed, 
and where it is properly treated, allowing one to assess how 
important this effect can be. Moreover, simulations assume 
sometimes choices for the IMF that are different than those 
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Figure 9. Comparison between different estimators of the stellar 
mass of galaxies, with respect to the real stellar mass, calculated 
as the simple sum of the mass of star particles Dashed 

and dot-dashed lines indicate simulations with a proper account¬ 
ing of stellar mass loss, while dotted lines are for simulations 
where mass loss is not consistently followed. 



used in the SPS models. It is for these reasons that we calcu¬ 
lated the stellar mass of simulated galaxies using the BC03 
method additionally to the direct result of the simulation 
(SIM method). 

In Fig. [ 9 ] we compare the stellar masses of simulated 
galaxies, obtained with the different methods, to the stel- 
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lar mass obtained directly from the simulations (SIM). Also 
shown is the 1-to-l relation (solid black line). In the lower 
panel, dotted, dashed and dot-dashed lines indicate respec¬ 
tively galaxies in the CS/CS"'"/MA samples, that differ in 
the treatment of mass loss of stars. We find that for galaxies 
where mass loss was not properly treated (CS), the stellar 
masses obtained with all methods are systematically lower 
compared to the direct result of the simulations. Further¬ 
more, the offset is very similar for all galaxies, of the order 
of 0.3 dex in log(M*). The results of the PHOTO, PETRO 
and MODEL methods are similar, with the former giving sys¬ 
tematically higher stellar masses, and the latter predicting 
the lowest values of stellar masses. 

In contrast, for galaxies in the simulations with a proper 
mass loss treatment (CS^/MA), the stellar masses obtained 
with the BC03 shows better agreement with the direct result 
(SIM), with differences up to ~ 0.1 dex. The SED method 
predicts slightly lower stellar masses, with typical differences 
of the order of 0.05 — 0.2 dex. In the case of the PHOTO, 
PETRO and MODEL methods, stellar masses have some scat¬ 
ter up to ~ 0.1 — 0.2 dex in log(M*), and are systematically 
lower for CS”*" (dashed lines) and higher for MA (dot-dashed 
lines) compared to the direct result of simulations. 

The results of this section show that the observational 
estimators are able to recover the mass up to 0.1-0.2 dex in 
logarithmic scale, only in simulations where the mass loss of 
stars is consistently included. In contrast, simulations with 
no treatment of mass loss can still be compared with obser¬ 
vations in a meaningful manner after a fixed correction (at 
least for the mass range of our sample) is applied. 


4.3 Stellar ages and metallicities 


The determination of stellar metallicities and ages of galax¬ 
ies from observational data is usually accomplished by fitting 
either selected sensitive absorption line indice s (e.g. Lick in¬ 
dices, |^^rthexet|aI][T99l, [Wage7^FaI][liil) derived from 
high S /N spectra to a nid of S PS rn odels fsee lTrager et al.l 
l2000alla iGallazzi et all l2005l . l2006li or the full available 
spectrum (A-to-A method, e.g. ICid Fe rnandes et al.l l2005l : 
iT^oieiro et al.l l2007l . l2009l : IChen et al.l 2010ll , which also re¬ 
quires high quality and S/N spectra for a reliable fit. Some 
authors also exploit colors to estimate ages and metal¬ 
licities, especially at high red shi ft (e. g. iLee et ah 1 I 2 OIOI: 
IPforr. Maraston fc Toniid l2012l : Ol2013f ). 

In this Section, we compare the stellar ages and metal¬ 
licities of the simulated galaxies obtained in different ways. 
In order to assess the effects of the fiber bias, for all methods 
we consider stars within the same field of view than in the 
previous section, i.e. in a projected area of (60kpc)^, and 
also only those within the fiber. The different methods are 
as follows: 


• SIM: the direct result of the simulations, i.e., the mean 
stellar age/metallicity of a galaxy is calculated as the cor¬ 
responding mass-weighted mean over the stellar particles. 
This requires no post-processing or additional assumptions, 
and is the most common way to derive ages/metallicities 
from simulations. 

• SIM-LUM: we weight with the luminosity of stellar 
particles to calculate the average ages/metallicities, which 
more closely reflects what would be obtained in an observa¬ 


tion; we use the r-band particle luminosity for stellar ages 
and the particles’ luminosity in all visible bands for the 
metallicities (both obtained with the BC 03 model) to mimic 
the SDSS analysis (IGallazzi et al.lBoOSh . 

• SED-FIT: we fit the spectra of simulated galax- 
ies obtained with sunris e using the starlight code 
dCid Fernandes et al.l I 2 OO 5 I I , selecting the best mixture of 
~ 300 instantaneous-burst SSPs with different ages and 
metallicities. The fitted spectra correspond to the face- 
on views and include dust and nebular emission (SRfaceon 
method). From the fitted mixture of SSPs we calculate the 
mean (optical) light-weighted ages and metallicities. 

• LICK-IND: the mean ages and metallicities were com¬ 
puted by Anna Gallazzi (priv. comm.) for a subsample of 
10/15 galaxies using the method described in IGallazzi et al.l 
(l2005h . The method is based on fitting five absorption 
features (Lick indices D4000n, H/I, H(5 a+H7a, [Mg 2 Fe], 
[MgFe]) to libraries constructed with the BG03 model with 
different star formation history, metallicity and velocity dis¬ 
persion, deriving the (optical luminosity-weighted) stellar 
metallicities and (r-band weighted) mean ages. As the spec¬ 
tra generated with SUNRISE have the resolution of the in¬ 
put stellar model (SB99) with spacing 20 A (which is too 
low to compute meaningful absorption indices), we thus ex¬ 
plore an alternative way to measure the Lick indices of the 
simulated galaxies directly from the BC03 tables, after in¬ 
terpolating with age and metallicity for each star particle 
and averaging weighting with the particle mass. However, 
we note that these may differ from the indices measured 
directly from the spectrum, as done in observations and in 
the model library of Gallazzi et al. used to interpret obser¬ 
vations. For this reason, this method is not fully consistent 
with the one used in SDSS. We will explore this point in 
subsequent analysis (Paperll) with higher-resolution spec¬ 
tra. 


The upper panel of Fig. [TO] shows the galaxies’ stellar 
ages obtained with the various methods described above as 
a function of the age obtained directly from the simulations 
(i.e. with the SIM method, Age®^”^), together with the 1:1 
relation, while the lower panel shows the corresponding dif¬ 
ferences. Results are shown in the cases of considering all 
stars in the field of view and stars within the fiber. 

In general, when all stars in the field of view are con¬ 
sidered (filled symbols and solid lines), the ages obtained 
using the SIM-LUM and SED-FIT methods are similar, and 
are systematically lower than the simple average of the age 
of stellar particles. This results from the fact that young 
stellar populations emit more light than old ones, therefore 
have a comparatively higher weight in SIM-LUM and SED- 
FIT. Typical variations are of the order of^ 2 Gyr, but the 
discrepancies are larger for some systems, with a maximum 
of about ~ 5 Gyr for the youngest galaxy. The estimation 
with the (mass-weighted) Lick indices gives systematically 
older ages for young galaxies (Age®^”^ < 8 Gyr) compared to 
SIM-LUM and SED-FIT, and slightly older ages compared 
to SIM, while it is in the range of the other observational 
estimators for the older galaxies. 

Gontrary to what we found for the ages estimated from 
all stars, in the case of the fiber quantities (open symbols 
and dotted lines) we find systematically higher ages com¬ 
pare to Age®^"^, with differences up to 4 Gyr. Note that 
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Figure 10. Difference between the mean mass-weighted stellar 
age in simulations and the ages calculated with the observational 
methods. The galaxies are sorted by increasing mean age to the 
right. 
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Figure 11. Difference between the mean mass-weighted stellar 
metallicity and the observational values, ordering the galaxies by 
increasing stellar metallicity to the right. 


the differences depend sensitively on the age but also on the 
presence/absence of age gradients which, as shown in Fig.[3l 
can vary significantly from galaxy to galaxy. For old galax¬ 
ies, all methods, in general, agree better with each other. 
We note that, for qnantities within the fiber, the two sources 
of bias (preferential sampling of inner populations and de¬ 
pendence on age gradients) can lead to both positive and 
negative differences with respect to the direct result of the 
simulation. The LICK-IND method taking only particles in 
the fiber (LICK-IND hber) gives in general slightly older 
ages compared to LICK-IND. 

We also applied the three methods and calculated total 
and fiber quantities for the stellar metallicities. Fig. llll shows 
a comparison of results. Compared to the direct result of 
the simulation, the SIM-LUM and SED-FIT methods give in 
general lower metallicities for metal-poor (old) galaxies and 
higher metallicities for more metal-rich (younger) ones. This 
is again explained by the different relative weight of old and 
young stars to the average metallicity (note that the spread 
in stellar metallicities is much smaller than that for stellar 
ages). Differences are however moderate, always of the order 
of ±0.3 dex. 

For quantities within the fiber, we also find moderate 
differences; in this case, metallicities tend to be higher com¬ 
pared to total quantities, as the contribution of very low 
metallicity stellar populations gets smaller (Fig. |4|). When 
we apply the LICK-IND and LICK-IND fiber methods, we 
hnd, for the low-metal sample {\og,{Z/Z q) < —0.45), system¬ 
atically higher metallicities, while for metal-rich galaxies we 


find better agreement with the other indicators. We detect 
small differences between LICK-IND and LICK-IND fiber, 
the latter giving in general slightly higher metallicities. It is 
notable that for the metal-poor galaxies, most observational 
methods disagree with the direct result of the simulations 
by an approximately constant factor, instead of scaling with 
the metallicity. 

Our results show that in our simulations it is not pos¬ 
sible to estimate the mean stellar age with accuracy bet¬ 
ter than ~ 3 Gyr, and the stellar metallicities better than 
~ 0.3 dex, depending on the method and on the way aver¬ 
aged quantities are calculated. The fiber bias has a strong 
effect on the derived ages/metallicities, due to the prefer¬ 
ential sampling of old and metal-rich stars; our 15 galaxies 
show wide variety of gradients, and therefore we find that 
the hber bias when applying different methods can signih- 
cantly change from galaxy to galaxy. 

4.4 Gas metallicity 

In this section, we compare the gas oxygen abundance of 
our simulated galaxies, using different estimators. Observa- 
tionally, the determination of the chemical composition of 
the gas in galaxies is based on metallicity-sensitive emission 
line ratios. A number of different metallicity calibrations are 
used to estimate the O/H ratio in nebulae and galaxies, and 
comparisons among the results of various calibrations re¬ 
veal large discrepancy wit h systematic offsets that can reach 
~ 0.7 dex (|Pilvuginll200ll . see also Appendix 

As in the previous sections, we calculate the gas chem- 
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ical composition following the most common methods from 
simulations and observations, and, when appropriate, con¬ 
sidering all gas particles and only those within the fiber 
(labelled as “hber”). When we use (face-on) sunrise spec¬ 
tra for the calculation (T04, KK04, Te), we first Balmer- 
correct the e mission line ratios for dust extinction using th e 
Calzetti law llCalzetti, Kinney fc Storchi-BergmannI Il994l . 
The methods are described in the following. 


• SIM: the gas metallicity of a galaxy is calculated as the 
mean (O/H) of gas particles, weighted by their mass. 

• HII: is the mean mass-weighted O/H of gas particles, 
but only considering particles around stars younger than 20 
Myr and inside 1 kpc radius. In this way, we mimic the 
preferential bias to young stellar populations in nebulae- 
based measurements. Note that although the HII regions 
have sizes which can vary in a wide range of 0.1 — 200 pc 
from ultra-compact to giant extragalactic HII regions, we 
are limited to define smaller sizes of the HII regions by the 
resolution of the simulations, which is of the order of ~ 1 
kpc. 


• T04 JTremonti et al.l|2004h : the gas metallicity is com- 
pu ted using the ca l ibrati on of the i? 23 -u pper branch given 


Tremonti et al.l (|2004l ~l. According to iKewlev fc EllisonI 


ll2008tl we use [Nil]/[OH] to remove the degeneracy in the 
i? 23 -metallicity relation, and we define the upper branch as 
Log ([Nil]/[ Oil]) ^ —1.2 (see also A ppendix lAll . 

• KK04 llKobulnickv fc KewlevI [2004l l : is a widely used 
metallicity calibration, based on an iterative method that 
allows to solve both the oxygen abundance and ionization 
parameter using the O 32 and 7?23 line ratios. 

• Te: the electron-temperature calibration (also referred 
to as “direct” method) is based on the ratio between the 
auroral line [OHI] A4363 and [OHI] A4959, 5007; it is com¬ 
monly used to determine gas metallicities when the weak 
auroral lin e can be detected. Here we follow the procedure 
outlined in lizotov et al.l (l2006l L 


Fig. [ 12 ] shows our results for the oxygen abundances 
of the simulated galaxies obtained using the different meth¬ 
ods. We show results as a function of the real gas metallic¬ 
ity, i.e. 12+log(0/H)®^'^, the differences are also calculated 
with respect to this estimator. As expected, the metallici¬ 
ties predicted using the SIM-fiber and T04-fiber methods are 
systematically higher compared to SIM as the most metal- 
enriched regions in the bulge are preferentially sampled in 
these cases. This happens however only for the metal-poor 
galaxies, while for the metal-rich ones these estimators agree 
very well. Note that, in our simulations, the galaxies with 
low/moderate metal content (CS and CS”*" samples) are also 
the ones with stronger metallicity gradients, while galax¬ 
ies in the MA sample do not show significant metallicity 
gradients, masking a possible bias due to preferential sam¬ 
pling (see also Fig. The HII estimator, which samples 
the gas near young star particles, shows ±0.5 dex scatter 
around the SIM value, in general higher for the galaxies with 
12±log(0/H)^ 8.3 and lower for the more oxygen-rich ones. 

The gas metallicities obtained with the emission-line 
calibration of T04 are systematically higher, with large dis¬ 
crepancies only for galaxies with low metal content. How¬ 
ever, note that only for 11/15 galaxies this calculation is 
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Figure 12. Comparison between the different estimators for the 
the galaxies’ oxygen gas metallicities, as a function of that ob¬ 
tained with the SIM method. For reference, in the upper panel 
we include the 1:1 relation as a solid line. As explained in the 
text, T04 calibration is valid only for a subsample of the galaxies. 


possibl43. Also for the KK04 method, we find good agree¬ 
ment with SIM for the metal-rich systems, with some offsets 
for metal-poor galaxies (i.e. 12 ± log(0/H) < 8.3). 

The largest discrepancies in gas metallicities, both with 
respect to the real value and to the other observational esti¬ 
mators, are found for the Te method, which gives the lowest 
metallicity values. Unlike the other methods, the systemat¬ 
ically lower metallicities are found even for the most metal 
rich galaxies where the rest of the methods agree well. The 
largest differences are of the order of ~ 1 dex, and only a 
weak trend with the metallicity is detected. 

It is worth noting that the emission line intensities in 
the spectra calculated from sunrise rely on the mappings 
HI photoionization code. The line ratios are then affected by 
uncertainties, assumptions and approximations in the model 
that describe th e photodissociation re gions, e.g. on-the-spot 
approximation dStasiriskal l2002l . l2007l f. assumptions on the 

J eometry of the nebula, dust grain c omposition, shocks, etc. 
Groves. Dopita fc Sutherlandll2003l . which will in turn af¬ 
fect the derivation of the gas metallicities. In any case, our 
results show that wide differences in metallicity can appear 
due to the use of various calibrations even on the same spec¬ 
trum modeled with a given photoionization code. 

In summary, our analysis shows that the fiber bias 
has a strong effect on the mean gas metallicity of galax- 


® Note that, when the T04 calibration is used inside the fiber, 
the estimation is possible for 12 out of the 15 galaxies. 
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ies, particularly for metal-poor systems and for systems with 
strong metallicity gradients, with metallieities systematically 
higher, up to ^ 1 dex, compared to the total mean metallic- 
ities. Deriving metallieities from the emission-line calibra¬ 
tions shows also large offsets, that can reach 0.5 — 0.8 dex 
compared to the value in the simulations, but in this case, 
the results are more diverse, with both positive and negative 
differences. Finally, the calibration based on electron temper¬ 
ature predicts systematically lower metallieities, by 0.7 — 1 
dex, compared to the rest of the methods and to the metal- 
licities derived directly from the simulations. 


4.5 Star formation rate 


In observations, the SFR of galaxies is estimated using dif¬ 
ferent methods at different redshifts, e.g. the luminosity of 
the Ha line, of the [0//]A3727 line and the UV continuum 
for low, intermediate and high redshifts, respectively. Each 
SFR indicator is affected by theoretical biases and obser¬ 
vational uncertainties (e.g. different timescale, sparse wave¬ 
length sampling, contribution from old stars, AGN contam¬ 
inatio n) which can influence the final SFR value llCalzettil 
I 2 OO 8 I ). Furthermore, the use of different SFR proxies at dif¬ 
ferent redshifts may introduce systematic errors that can 
bias the comparison between simulations an d observations 
llRosa-Gonzalez. Terlevich fc Terlevichll2002l ). 

In this section, we calculate the SFRs of our simulated 
galaxies with differ e nt m ethods. We use the calibrations 
given in iKenniciifrO (Il998l) , app lying the correction factor 
/imp = 1.5 ( Calzetti et al. 200^) to account for the use of 
a different IMF ( Kennicutd Il998l assumes a Salpeter IMF 
while we assume a Kroupa IMF in sunrise and Ghabrier 
in BC03). We note that, in SDSS, the fiber-derived SFRs 
llBrinchrnann_et ah 20041 ) are corrected to total SFR with 
the Salim et al. ll2007l ) method. For this reason, we use the 
full face-on spectra obtained with sunrise, dust-corrected 
with the Galzetti law, to derive SFRs, instead of the spectra 
within the fiber. The various methods to derive SFRs are 
described below. 


• SIM: the real SFR extracted directly from the simula¬ 
tions (stellar mass formed over time interval), averaged over 
the past 0.2 Gyr. 

• BC03: we convert the rate of ionizing photons calcu¬ 
lated with BC03 into SFR: 

SFR(Moyr-i) = 1.08 /jmV x 10-®®Q(77°) (s’^). 

• Ha: we extract the Ha-luminosity from the sunrise 
spectrum and convert into SFR according to: 

SFR(Mq yr-i) = 7.9 /imV x (erg s'^). 

• UV: we calculate the flux from sunrise spectra of the 
nearly-fiat region between 1500 — 2800 A and calculate the 
SFR as: 

SFR(M 0 yr-') = 1.4, ff^^ x lO""®!,. (erg s'^Hz-^). 

. [Oil]: we estimate the SFR using the luminosity of the 
forbidden-line doublet [OHJA3726, 3729 in SUNRISE: 

SFR(Mo yr-i) = 1.4 f{^^ x 10-"L(OH) (erg s'^). 

• FIR: the SFR is estimated from the sunrise FIR lumi¬ 
nosity integrated over the range 8 — 1000 /rm as: 

SFR(M 0 yr-i) = 4.5 /jmf x IQ-^^Lfir (erg s'^). 
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Figure 13. Comparison between the SFRs obtained using var¬ 
ious methods as a function of the real SFR of the simulation 
(method SIM). In the upper box, we also include the 1:1 relation 
for reference. 


All observational methods are sensitive to the emission 
from young massive stars, although the UV indicator has 
10 times longer timescale (~ 100 Myr), because massive 
stars stay luminous for longer time in the U V with small 
production of ionizing photons llGalzettill2008l) . 

Fig.[T3]shows the differences between the SFR obtained 
with the different estimators and that obtained directly from 
the simulations (SFR®^^). Except for the galaxy with the 
lowest SER, the SFRs obtained with BC03, Hq, and UV 
methods agree well, with differences of < 0.2 — 0.4 dex in 
logarithmic scale with respect to the real SFR of the sim¬ 
ulations. Note that BC03 and H^ are based on the same 
method (conversion of the number of ionizing photons into 
SFR), but use two different models. 

The other observational indicators ([OH] and FIR) pre¬ 
dict in general lower SFRs compared to SFR®^*^, with max¬ 
imum differences of about 1 dex. In particular, the [OII]- 
SFRs of metal-poor galaxies are significantly lower than 
SFR®^^, by ~ 0.5 dex. Note that the [Oil] line inten¬ 
sity and [Oil]/Ha ratio strongly depend on the metallic¬ 
ity ( Kewlev. Geller fc Jansenll2004h and has been calibrated 
byi ennicu 13 lll998l) for around-solar metallicity; this cal¬ 
ibration is not accurate for the metal-poor galaxies of our 
sample. Also in the case of the SFR obta ined with t h e FIR 
method differences are expected, as the iKennicii tl (Il998l) 
calibration is valid only for dusty metal-rich starburst galax¬ 
ies (note that calibrations that follow both dust-obscured 
and uno bscured star f ormation can be found in the litera¬ 
ture, see lGalzettill2008h : for the metal-poor sample, since the 
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amount of dust in sunrise is proportional to the metallicity, 
we have only small dust absorption ISec. 14.41) and reemission 
in the infrared (see also lHavward et alfeOld l. To further in¬ 
vestigate this effect, we plot in Fig. [13] the result obtained 
using the FIR method of one of our metal-poor galaxies (C- 
CS, orange square) where the assumed dust-to-metal ratio 
in SUNRISE was increased 10 times (as in Section[4T}. In this 
case, the SFR raises by ~ 0.6 dex, although it is still moder¬ 
ately lower compared to the results of the other estimators. 

Finally, note that SFRs derived observationally and in 
simulations assume different timescales, of about 10 — 100 
Myr for the former compared to 0.2 — 0.5 Gyr for the latter. 
This can result in signihcant discrepancies among the two 
estimations if we are in the presence of recent starbursts. In 
our sample, however, we observe such a systematic difference 
only in the galaxy with the lowest SFR (D-CS). 

The results of this section show that the SFRs indica¬ 
tors exhibit differences of < 0.4 dex in log(SFR), which are 
caused both by the difference in methods and the different 
characteristic time-scales to which the methods are sensitive. 
The Ha and UV methods predict similar SFRs compared to 
the direct results of the simulations, except for the galaxy 
with the lowest SFR. For galaxies with low metallicities the 
[on] method gives systematically lower SFRs, with absolute 
differences lower than 0.5 in logarithmic scale. For the FIR 
method the differences are larger in metal-poor galaxies, with 
SFRs systematically lower by 0.5 — 1 dex in log(SFR), while 
the FIR estimation agrees with the SFR of the simulations 
up to 0.1 — 0.2 dex for metal-rich galaxies. 


5 DISCUSSION AND CONCLUSIONS 

We have used a set of 15 simulated galaxies, of similar 
mass to the Milky Way, to study biases and systematics 
in the derivation of observables, and in the comparison be¬ 
tween them and observational data. The aim of this study is 
twofold; first, help simulators to be aware of systematics and 
be able to reliably judge the agreement between simulated 
and observed galaxies and, second, help observers to inter¬ 
pret observational results by being able to better quantify 
the differences between the observationally-obtained galaxy 
properties and the real ones (that are known in the simula¬ 
tions). 

Our simulations comprise 15 galaxies with a variety of 
merger, formation and accretion histories, which results in a 
variety of final morphologies, star formation rates, gas frac¬ 
tions and metallicities. As these properties somewhat de¬ 
pend on the modelling of feedback, the same 5 galaxies were 
simulated using three different models for chemical/energy 
feedback; the strength of feedback and the amount of chem¬ 
ical yields varies from moderate to strong. For moderate 
feedback we find galaxies that are more metal-poor and form 
their stars earlier, while for strong feedback we find younger 
galaxies with a higher metal content. This diversity, typical 
of real galaxies, makes the sample ideal to test biases and 
systematics in the derivation of the synthetic spectra, and 
to study the dependence of such effects on galaxy properties 
(with the caveat that with our sample we do not represent 
the more metal rich galaxies observed). 

For our study we have computed the synthetic spectra 
of the simulated galaxies at z = 0, as at low-redshift many 


databases of galaxy properties derived from large galaxy 
surveys (e.g. SDSS, 2dFGRS and 6dFGS) are available. 
For this purpose, we followed three different approaches: 
(i) Stellar Population Synthesis (SPS) models, which give 
the spectra coming from stars; (ii) SPS models including 
dust extinction with a simple recipe; and {Hi) a full ra¬ 
diative transfer calculation that gives the spectra including 
stellar and nebular emission, as well as the effects of dust 
which are parametrized using the metallicity of the interstel¬ 
lar medium. We have used the synthetic spectra to derive 
the observables (magnitudes/colors, stellar masses and ages, 
stellar/gas metallicities and star formation rates) in various 
ways, as to mimic real observations, and we have compared 
the results with the direct outputs of the simulations, i.e. 
the real galaxy properties. 

Biases and systematics appear at various stages in the 
process of obtaining the observables from the simulations, 
due to: 

- assumptions and parameters of SPS models, 

- dust/radiative transfer/projection effects, 

- weighting with the mass instead of luminosity for mean 
quantities, 

- observational biases, such as extrapolation to external re¬ 
gions of galaxies where no spectral/photometric data are 
available (e.g. fiber bias, Petrosian/Model magnitudes), 

- different parametrization of the star formation history and 
dust extinction when quantities are derived fitting a pre¬ 
constructed grid of models, 

- in the particular case of gas metallicities and SFRs, the 
use of different calibrations. 

We tested the effects of such biases on the magnitudes 
and colors of simulated galaxies, their stellar masses, ages, 
stellar/gas metallicities and star formation rates. Our results 
can be summarized as follows: 

• Magnitudes 

- The galaxies’ magnitudes in the {u,g,r,i,z) SDSS bands 
derived from the five different SPS models used in our work 
show good agreement, with differences lower than 0.1 dex. 
Differences are larger for the most metal-rich galaxies, which 
also have a higher contribution of young and intermediate 
age stars, for which model uncertainties are the largest. 

- Dust effects can be important if galaxies are seen edge-on. 
A simple angle-averaged dust model predicts galaxies that 
are ~ 0.4 — 0.8 dex fainter compared to the no-dust case. If 
full radiative transfer is considered, where dust is traced by 
the metals, edge-on galaxies appear in general 0.3 — 1 dex 
fainter, but only if these are relatively metal-rich. 

- Estimating the magnitudes using a more observational ap¬ 
proach (as for instance the Petrosian and Model magnitudes 
of SDSS) exhibits some differences compared to the real 
magnitudes of the simulated galaxies, with offsets ~ 0.2 —0.3 
dex for 60% of the systems and up to 0.6 dex for the remain¬ 
ing 40%. 

• Stellar masses: 

When observational techniques are applied to the simu¬ 
lated galaxies to estimate their stellar masses, the results 
vary depending on the treatment of mass loss in the simu¬ 
lations: 

- If mass loss includes only SNe (which is a small effect), 
all observational estimators give systematically lower stellar 
masses compared to the direct result of the simulation, as 
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the fitted models include the full mass loss of Stellar Popu¬ 
lations at each stage of evolution. The offset is similar for all 
galaxies, of about a factor 0.3 in log(M*) or, equivalently, of 
50% in M*. 

- If mass loss is properly treated in the simulations (e.g. 
adding AGB stars), the observational techniques recover the 
real stellar masses with differences 0.1 — 0.2 dex in log(M*). 

• Stellar ages: 

- If stellar ages are estimated mimicking observational tech¬ 
niques weighting with the luminosity (but ignoring the fiber 
bias), galaxies appear younger, typically by ~ 2 Gyr, com¬ 
pared to the direct result of simulations (mean age of stellar 
particles). Some galaxies exhibit larger differences, up to a 
maximum of 5 Gyr; among them, we find both very young 
and very old galaxies. 

- If only stars in the nuclear part of galaxies are in¬ 
cluded (as to mimic single-fiber surveys such as SDSS), the 
observationally-estimated galaxy ages are in general higher 
compared to the direct result of the simulations with dif¬ 
ferences of about 2 — 4 Gyr for young systems and of < 1 
Gyr for old ones. However, the results depend sensitively on 
the particular properties of galaxies, namely the presence 
of significant/insignificant age gradients and the resulting 
preferential sampling of old stars within the fiber. 

• Stellar metallicities: 

- Observationally-estimated stellar metallicities (ignoring 
the fiber bias) are in general lower compared to the di¬ 
rect result of the simulations, in particular if galaxies are 
metal-poor, with differences up to —0.3 in logarithmic scale. 
For more metal-rich systems, differences are similar in ab¬ 
solute value, but can vary from positive to negative. The 
differences in the behaviour of metal-rich/metal-poor galax¬ 
ies originates in the different weight of old/metal-poor and 
young/metal-rich stars in the different methods, in particu¬ 
lar depending on the way the mean metallicity is calculated 
(mass/luminosity-weighted). 

- If only stars within the fiber are considered, stellar metal¬ 
licities tend to be higher compared to the direct result of the 
simulations with typical differences of ~ 0.1 — 0.3 dex and 
no strong dependence with the real stellar metallicity. 

• Gas metallicities: 

- Deriving gas metallicities using different emission-line cal¬ 
ibrations show large spread in the resulting metallicity val¬ 
ues, up to 0.5 — 0.8 dex. In general, metallicities lower than 
the real value are predicted for very metal-poor systems, 
while for more metal-rich galaxies they show a better agree¬ 
ment. Significant differences are found when the calibration 
based on electron temperature is used; in this case the de¬ 
rived metallicities are systematically lower by differences of 
the order of 0.7 — 1 dex. 

- Gas metallicities are systematically higher than the di¬ 
rect result of simulations when the fiber bias is included, 
due to a preferential sampling of metal-rich regions in the 
galaxy. Differences can be up to 1 dex for the most metal- 
poor galaxies, while the for metal-rich sample differences are 
always smaller than 0.1 dex. 

• Star formation rates: 

- Observationally-derived SFRs using different indicators 
present in general offsets of the order of ±0.4 dex in 
log(SFR), mainly due to the differences in methods and 
time-scales to which they are sensitive. The largest differ¬ 
ences are found for the estimations based on the FIR, which 


gives systematically lower SFR values up to 1 dex in loga¬ 
rithmic scale for the metal-poor galaxies. 

In summary, we have shown that it is important to prop¬ 
erly take into account different observational biases when 
galaxy observations are interpreted and linked to different 
underlying physical processes. Furthermore, a meaningful 
comparison between observations and simulations of galax¬ 
ies requires a good understanding of the systematics; if this 
is not considered it is not possible to properly judge agree¬ 
ment/disagreement between simulations and observations, 
and to decide which of the physical processes included in 
the simulations are more relevant in the context of the for¬ 
mation and evolution of galaxies in a cosmological context. 

In a companion paper (Guidi et ah, in prep.), we com¬ 
pare the properties of our simulated galaxies to observa¬ 
tions of the SDSS survey, in order to gain more insight in 
their agreement/disagreement in relation to the treatment 
of feedback, star formation, mass loss and chemical enrich¬ 
ment in the simulations. We also quantity differences from 
the observationally-derived quantities and the direct result 
of the simulations to provide the relevant scalings and error 
bars for a meaningful comparison, in the case of the different 
galaxy properties. 

Finally, we note that other possible sources of biases and 
systematics have not been explored here, e.g. related to the 
inability to resolve the typical height of gaseous discs (which 
could affect the comparison between face-on and edge-on 
projections). Moreover, how well the simulated galaxies re¬ 
produce the observed sizes/concentrations/luminosity pro¬ 
files can affect our findings related to the fiber bias and 
Petrosian/Model Magnitudes (although in previous papers 
we have shown that our models produc e galaxies in broad 
agreement with observational results, e.g .lScann anieco et al.l 
l2010l . IScannapieco et aLll2012l . lAumer et al.ll20l3h We hope 
that our work encourages other researches to explore all 
these issues, as well as we will continue to test systemat- 
ics/biases in further analysis. 
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APPENDIX A: GAS METALLICITY 
CALIBRATIONS 


The observational derivation of the gas metallicity of galax¬ 
ies requires a calibration of emission lines to recover 
the (O/H) ratio; different ones are used in various ob¬ 
serva ti onal campaigns (e.g Z^ritgkyjKgnnicutI_^H^£hr^ 
]J19^;_ jKobulnlc' KewlevI 20041 : Kewlev et ah 20021 : 
Pettini & Pagel 2004h . Previous work has already shown 


that large discrepancies, as large as ^ 0.7 dex, arise when 
different calibrations are used dPilvuginl [20o3). The cal¬ 
ibrations can be broadly classified as empirical (further 
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Table Al. Main characteristics of the different metallicity cali- ♦ A-CS^ ^ B-CS^ ■ C-CS^ • D-CS^ ^ E-CS^ 

brations fsee lKewlev fc EllisonllioO^ for more details'). A-CS I B-CS M C-CS | D-CS b E-CS 

+ A-MA • B-MA w C-MA 4 D-MA 4- E-MA 


Method Class Reference 


T04 

theoretical 

Tremonti et ah l'2n04') 

KK04 

theoretical 

Kobulnickv & Kewley 12004') 

Te 

direct 

Izotov et ah ([2006) 

Z94 

theoretical 

Zaritskv. Kennicutt Sz Huchra fl994) 

KD02 

theoretical 

Kewley et al. (2002) 

M91 

theoretical 

^cGaugh 119911 

D02 

combined 

Denicolo. Terlevich & Terlevich 120021 

PP04 

empirical 

Pettini & Pagel 120041 

PP04(O3N2) 

empirical 

Pettini fc Pagel 120041 

P05 

empirical 

Pilvugin. 'Vflchez & Thuan (20101 


subdivided into “direct” and “statistical”) and theoreti¬ 
cal. The former derives metal abundances directly from 
electron temperature-sensitive lines or relations between 
temperature-estimated metallicities and strong emission 
lines, while the latter relies on photoionization mode ls to cal¬ 
ibrate the metallicity ind icators (see the reviews bv iFerlandl 
I 2 OO 3 I : IStasihskal [20011 and iKewlev &: Ellisonll200^ '). 

In general, theoretical calibrations give, for a 
given galaxy, a higher metallicity compared to elec¬ 
tron temperatur e -based estimators llLiang et al.l l2006l : 
iKewlev fc EllisonI I 2 OO 8 II . The causes of these discrep¬ 
ancies are still unclear, although some authors con- 
.i ecture either problems in the ph otoionization models 
llKennicutt, Bresolin fc Garnett! I 2 OO 3 II. or temperatu r e gra¬ 
dient fluctuat ions in the nebulae (IStasihskaf l2005l . l2007l : 
lBresolir]|2008l l. The effects of these disagreements can also 
bias the determination of th e shape and the zero poin t in 
the mass-metallicity relation ([Andrews k, Martinill201^ . 

Fig. presents a comparison between the oxy¬ 

gen abundance of our simulated galaxies assuming var¬ 
ious calibrations, that we summarize in Table ET] All 
results are derived after Balmer-correcting the emission 
line ratios for dust extinction using the Calzetti law 
dCalzetti. Kinney fc Storchi-Bergmannlll994l l. Some calibra¬ 
tions have limited ranges of validity, and here we only 
show results for galaxi es that satisfy these conditions (see 
iKewlev fc EllisonI 120081 for details). 

It is clear that the metallicities for the various emission¬ 
line calibrations show larger differences, up to ~ ±1 dex. 
The use of different calibrations is certainly a source of con¬ 
cern, although in general the discrepancies are similar for 
all galaxies regardless their metallicities, although variations 
are smaller for the more metal-rich galaxies. In particular, 
the calibration based on electron temperature is in general 
lower than the others (see Sec. 14.411 . manifesting itself as a 
clear offset, particularly for galaxies with (direct) metallic¬ 
ities higher than 12 -|- log{0/H) ~ 8.2. On the other hand, 
T04 and Z94 give the highest metallicity values, with the 
other calibrations lying in between. For the most metal-rich 
galaxies, a number of indicators agree relatively well, par¬ 
ticularly KK04, KD02, M91, and even T04 and Z94, with a 
small spread of 0.1 — 0.2 dex (note that all of these are 
theoretical calibrations). For the most metal-poor galaxies, 
we also find that various of the calibrations agree relatively 
well with each other, still with variations of ±0.5 dex. 
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Figure Al. Comparison of the gas metallicity of simulated galax¬ 
ies using different calibrations, as a function of the direct result 
of the simulation. 
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